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Abstract. In the case of barotropic FRW cosmologies, the Hubble parameter in conformal 
time is the solution of a simple Riccati equation of constant coefficients. We consider these 
cosmologies in the framework of nonrelativistic supersymmetry that has been effective in 
the area of supersymmetric quantum mechanics. Recalling that Faraoni [Amer. J. Phys. 67 
(1999), 732-734] showed how to reduce the barotropic FRW system of differential equations 
to simple harmonic oscillator differential equations, we set the latter equations in the su- 
persymmetric approach and divide their solutions into two classes of 'bosonic' (nonsingular) 
and 'fcrmionic' (singular) cosmological zero-mode solutions. The fermionic equations can be 
considered as representing cosmologies of Stephani type, i.e., inhomogencous and curvature- 
changing in the conformal time. We next apply the so-called shifted Riccati procedure by 
introducing a constant additive parameter, denoted by S, in the common Riccati solution of 
these supersymmetric partner cosmologies. This leads to barotropic Stephani cosmologies 
with periodic singularities in their spatial curvature indices that we call U and V cosmolo- 
gies, the first being of bosonic type and the latter of fcrmionic type. We solve completely 
these cyclic singular cosmologies at the level of their zero modes showing that an acceptable 
shift parameter should be purely imaginary, which in turn introduces a parity-time (PT) 
property of the partner curvature indices. 
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1 Introduction 

Riccati nonlinear equations have many applications in physics and their solutions under the 
name of superpotentials play an important role in supersymmetric quantum mechanics pQ (for 
a chronological review see [2]). In the cosmological framework, the supersymmetric methods 
have been recently reviewed in the book of Moniz [3] and in general the occurrence of Riccati 
equations is well documented, see for example [3]. With the cosmological application that we 
will consider in this paper in mind, we write the Riccati equation in the form 

R' + cR 2 + f = 0, (1) 

*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Benasque, Spain). The full collection is available at 
|http://www.emis.de/journals/SIGMA/SUSYQM2010.html1 
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where c is a real constant and / is a function of the independent variable. It is well known 
that particular solutions of the Riccati equation enter the factorization brackets of second order 
linear differential equations (usually with c = 1) 

4- + cR ) f 4- ~ cR ) u = & u" - c(cR 2 + R')u = ^ u" + cfu = 0. (2) 
at J \dt J 

The connections i? = or u = e c ^ ^ between the particular solutions of the two equations 
are also basic results of the factorization method together with the construction of the so- 
called super symmetric partner equation of equation ([2]) obtained by reverting the order of the 
factorization brackets: 

(it ~ cR ) (jE + cR ) v = v " ~ c ^° r2 ~ R '^ V = v " + C ^ + 2R '^ V = °' 
2 Shifted Riccati procedure 

In previous works, we have introduced a supersymmetric procedure in which we examined the 
consequences of a constant shift of the Riccati solution [U [6], i.e., 

R s (t) = R(t) + S. 

Rs obeys a Riccati equation of the form: 

R' s - 2cSR s + cR 2 s + (/ + cS 2 ) = 0, 

which for S = turns into the simple form ([1]). The corresponding linear second-order differential 
equation obtained by substituting Rs = is: 

u" s - 2cSu' s + (cf + c 2 S 2 )u = 

and one can immediately see that a particular solution is us = e cSt u. Similarly, vs = e cSt v is 
a particular solution of the supersymmetric equation 

v's - 2cSv's + c(f + 2R' + cS 2 )v = 

related to the non-standard Riccati equation 

-R's - 2cSR s + cR 2 s + (/ + 2R' + cS 2 ) = 0. 

However, if we consider the shifted Riccati solution Rs directly in the factorization brackets 
we obtain the following pair of supersymmetric equations: 

+ cRs ) (jE ~ cRs ) u = ° ^ U " + C U ~ 2c SR ~ cS 2 )U = 0, (3) 

(jt ~ ° Rs ) (jt + cRs ) V = ° V " + c (f~ 2cSR ~ cS2 + 2R "> V = °- 

These equations reduce to the equations for u and v when S = 0. Nevertheless, their solutions 
are not connected to the solutions u and v in the same simple way that us and vs are. In 
the rest of the paper, we will present in full detail the solutions of the latter equations in the 
interesting cosmological case of barotropic FRW cosmologies. 
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3 Cosmological Riccati equation in FRW barotropic cosmology 

In 1999, Faraoni [7] showed that the well-known comoving Einstein-Friedmann dynamical equa- 
tions of barotropic FRW cosmologies can be turned into a single Riccati equation for the Hubble 
parameter in conformal time T~L(r]) of the form 

H' + lU 2 + «7 = 0, (5) 

where henceforth the prime and also J- stand for the derivative with respect to rj, 7 = |7 — 1 is 
related to the adiabatic index 7 of the cosmological fluid and k = ±1 is the curvature parameter 
for the closed and open universe, respectively. Notice that for 7 = | equation © reduces to 
Ti' = which is not a Riccati equation and therefore the arguments of this paper do not apply 
in this case. Moreover, equation ([5]) is just the Riccati equation of the harmonic oscillator: 

R + oj R 2 + w = 0, 

if one sets 7 = ujq for the closed universe case. For the open case the analogy is with the up-side 
down harmonic oscillator. The main difference resides merely in the independent variable, which 
in the cosmological framework is the conformal time whereas in the harmonic oscillator case is 
the usual Newtonian time. In the following, we will consider only the k = 1 case since it allows 
us to focus on the periodic features of the problem. 

Usually the conformal Hubble parameter is defined as the logarithmic derivative of the scale 
factor of the universe, T-i{rj) = — , which is related to the observational (comoving) Hubble 
parameter through %{rj) = a(t)H(a(t)), for a discussion of this issue see [SJ p. 96]. Here, we 
define the conformal Hubble parameter directly in terms of the zero modes, i.e., the particular 
solution u, as = Substituting % u in equation ([5]) leads to the very simple (harmonic 

oscillator) second order differential equation 

u" + k u 7 2 u = 0, k u = 1. (6) 

Employing a terminology used in supersymmetric quantum mechanics we have previously 
called the u modes as bosonic zero modes [6]. They are 7 powers of the scale factor parame- 
ters a(rj), i.e., u = cp(ri). Faraoni used the following particular solutions: 

u\ ~ cos7?7 ai u (r]) ~ u 1 ^ . 

Using %i\ in the definition of l~Lu 

one gets 7i Ul = — tan 777. 
A class of cosmologies with inverse scale factors with respect to the standard barotropic 
ones but with a specific conformal time dependence of the curvature index can be obtained by 
considering the supersymmetric partner equation of equation ([6]), which is given by [6] 

v" + n v (i])^ 2 v = 0, 

where 

^(77) = -(1 + 2 tan 2 777) (7) 

denotes the conformal time dependent supersymmetric partner curvature index of fermionic type 
associated through the mathematical scheme to the constant bosonic curvature index. 
A particular fermionic solution v is of the following type 
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We can see that the u and v barotropic cosmologies are dual to each other from the standpoint 
of these particular solutions, in the sense that u\V\ = 7 and therefore the geometric mean of 
their scale parameters 

a g = (a lu a lv ) 1/2 = f^) 1 / 2 ^ 

is constant. Thus, a joint evolution of a u cosmology of constant curvature index and a v 
cosmology of the time-dependent curvature index (JT]) is stationary in conformal time from the 
standpoint of their geometric mean scale parameter a g . Since the fermionic cosmology is of 
time-dependent (periodic) curvature index, it can be also called a Stephani type cosmology [9] 
(for a review see |10|). although in conformal time. The fermionic metric is of the form: 



ds =a lv {rj) 



- dr ] . = +rW 



1 - n v (rj)r^ 



This metric should be thought of as an averaged metric in an inhomogeneous cosmology and so 
it does not even have to satisfy the Einstein field equations [TT] , although we prefer to keep this 
constraint for philosophical reasons (that is why we consider to be of Stephany type). The only 
thing it has in common with the standard FRW model is the Riccati solution. Such metrics, 
with other time-dependent curvature indices have been used to mimic the backreaction of small 
scale density perturbations on the large scale spacetime geometry 112]. 



4 Barotropic conformal-time Stephani-type cosmologies 
from a constant shift of the Hubble parameter 

We would like now to investigate the implications of the constant shifting S of the conformal 
Hubble parameter, according to the procedure discussed in Section [2j This means that we will 
build Schrodinger equations based on the Riccati solution Hs(v) = Huiiv) + S, similar to the 
equations ([3|) and ((4|) 

drj J \dr] J 

U" + K u {T])fU = 0, Kuirj) = 1 -5 2 + 25tan7? ? (8) 

and 

T- ~ ^ Us ) (-T + ^ ns ) V = ^ 
dr/ J \dr] J 

V" + lC u {r])fV = 0, M7 ? ) = -l-,S 2 + 2Stan77 ? -2tan 2 7? ? . (9) 

Equations ([8]) and ([9]) define two new classes of barotropic-like cosmological universes that 
we call hi and V cosmologies, respectively. One can also write averaged conformal-like metrics 
with the corresponding variable curvature indices. These two modified cosmologies are periodic, 
of period T = 5, and have the same conformal Hubble parameter given by: 

H u (v) = Hviv) = n s (v) = + s. 

cos 77/ 

Obviously, "Kuy —> W«„ when S — > 0. Both cosmologies have periodic singularities in their 
time-dependent curvature indices. In addition, Hs(0) = S, i.e., S serves as an initial condition 
for the conformal Hubble parameter of these cosmologies. 
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The linear independent solutions U\ and U2 have the following form 

U^n) = e~ iQ ^ 2 iq (1, -iS; 2 - iS; - e " 2 ^) , (10) 

or Wi(r?) = e~ s ^ [z 2 F 1 (l, -iS; 2 - iS; -z 2 )] , where z = e~^, Q s = lO—iS) and 2 F\ (a, b; c; z) 
is the hyper geometric function, and 

U 2 (rj) = e^ (2 cos 2 777-?; sin 27?/) ~e 5 ^ cos 777. (11) 

We also notice that the simple convergence condition ft(a + b — c) < for the hypergeometric 
series is fulfilled for all real values of 77. 

On the other hand, the linear independent solutions Vi and V 2 are given by: 

Vi(r?) = -„ (12) 

(2 cos z 777 — 7 sin 2777) 

or = #^ and 

V 2 (t/) =7 2 e ^s'' [(l-itan7r ? ) + i5(l-i5)(2cos 2 7r ? -isin27r ? ) - 2iS] , (13) 

or V2O7) = 27 2 e 5 ^ 2 cos 7?? + ^(^ cos 7^? + sin iv) ■ 

One can notice immediately that we have written the solutions (|10p . (jlip and (|12p . (|13p in 
a convenient Floquet-Bloch form that allows us to make the following comments. The parame- 
ter S affects only the period of the phases ^ of the solutions but not that of their periodic 
part. The solutions (JTOJ), and the Bloch factors e ±iUsri in (JT2J), (USD are bounded if and only 
if the "quasifrequency" has a real value, or equivalently 

j(l -iS) E ft. 

Taking into account that 7 € ft, the last condition implies that S is pure imaginary, i.e., 5 = is, 
s G ft. Thus, equation([8]) has bounded solutions Vs G ft. Notice that for a purely imaginary 
shift parameter, the curvature indices fC u and K, v are related through JC u>l/ (—rj) = /C* u (r/), where 
* denotes the complex conjugation operation. This means that we can have in this cosmological 
context the parity-conformal time (PT) symmetry. Additionally, by inspecting the solutions 
(flop . (fTTI) and (fl"2j) . (fT3j) we note that they are periodic for s p = (2m — 1) , m = 0, 1, 2, . . . and 
antiperiodic for s a = 2m , m = ±1, ±2, .... Plots of the solutions in the periodic case m = 1 are 
displayed in Figs. Q] and [2j whereas the antiperiodic case m = 2 can be found in Figs. [3] and fU 
In general, the IA modes are regular indicating that these shifted cosmologies are not sensitive to 
the singularities of their curvature indices at the level of their zero modes. On the other hand, 
the V cosmologies have periodic singularities in their imaginary parts but not in their real parts. 
In addition, the duality property is maintained for the pair of zero-modes U 2 V\ = const. 
The periodic case for m = gives the solutions 

Ul (77) = 2%e*™ sin 777, (77) = 2e~*™ cos 777 

and 

Vffa) = V§fo) = 1^- - 2 7 2 - 

2 cos 77? cos 777 

The solution U\{rj) is further discussed in the Appendix. 

Regarding the dualities we introduced here, it is worth emphasizing certain similarities with 
the superstring dualities [13] and the phantom duality |14| . In the first case, there is an invariance 
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Figure 1. The real (solid lines) and imaginary (dotted lines) parts of the periodic zero modes Ui(r]) (in 
red) and U2(j]) (in blue) for the shift parameter S = 3i. 
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Figure 2. The real (solid line) and imaginary (dotted lines) parts of the periodic zero modes Vi(r]) (in 
red) and V2{r]) (in blue) for the shift parameter S = 3i. 
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Figure 3. The real (solid lines) and imaginary (dotted lines) parts of the antipcriodic zero modes Ui(ri) 
(in red) and Ui(j]) (in blue) for the shift parameter S — Ai. 

of the action with respect to an inversion of the cosmological scale factor and special shifts of 
the value of the dilaton field (a a -1 and $ — > — 6 In a, respectively). We also have the 
scale factor duality and we also considered a shift but of the conformal Hubble parameter. 
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However, our context and mathematical techniques appear to be different. Moreover, the string 
action with such properties corresponds to cosmologies that are spatially flat and homogeneous. 
Thus, our scale factor duality connecting homogeneous and inhomogeneous nonflat cosmologies 
look more general. On the other hand, Dabrowski et al. discuss the symmetry 7 —¥ —7 of 
a nonlinear oscillator equation. In our paper, the 7 parameter occurs mostly trigonometrically. 
For the unshifted bosonic and fermionic cosmologies, we see readily that 7 — > —7 is a symmetry 
preserving their curvature indices, although we have now u±vi = —7 and a g = (— 7) -1 ^ 2 ^, 
a different constant. In the case of the shifted cosmologies, the sign changes 7 — > —7, K. u {rj; S) — > 
JC u (r]] —S), /C„(r?; S) — > /C„(r/; — S) leave the equations unchanged. 

5 Conclusion 

In this paper, we have introduced two classes of supersymmetric-partner classes of barotropic 
cosmologies of Stephani type, i.e., of variable spatial curvature index, by considering a constant 
shift of the conformal Hubble parameter of FRW barotropic cosmologies. This has been worked 
out in a super symmetric context similar to super symmetric quantum mechanics and can be 
applied to any type of cosmological fluid unless 7 = |. In other words, these new classes of 
inhomogeneous cosmologies together with the unshifted fermionic cosmology can be associated to 
any barotropic cosmological fluid with the exception of the coasting (non-accelerating and non- 
decelerating) universe. Interestingly, we have found that an acceptable shift is purely imaginary 
and moreover serves as an initial condition for the conformal Hubble parameter. (It is known 
that 'unnatural' initial conditions are required to explain why a thermodynamic arrow of time 
exists.) Cyclic behavior in conformal time of the cosmological zero modes can be obtained 
in addition to that of their curvature indices that in the pure imaginary case are also related 
through the parity-time (PT) property. 

On the other hand, if we forget about the curvature-changing universes, the same results 
can be interpreted as due to the time dependence of the adiabatic indices of the cosmological 
fluids |15j . Along this line, we notice that Dabrowski and Denkievicz [16] have provided a recent 
discussion of an explicit barotropic model in which the cosmological singularity occurs only 
in the singular time-dependent barotropic index. They assert that physical examples of such 
singularities appear in f(R), scalar field, and brane cosmologies, see [17] for a recent review in 
relation to all sorts of singularities. We think that the barotropic inhomogeneous cosmologies 
with periodic singularities in the curvature index that we introduced here could be related to the 
same physical examples. Moreover, the periodic singularities are not an impediment to build 
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appropriate cosmological zero modes along the whole conformal time axis that can be used 
to define novel classes of scale factors corresponding to these generalized barotropic universes 
within any chosen period of the curvature index. 

A Further remarks on the solution tl\{rf) 

Formula (15.3.4) in [18] (in Linear transformation formulas, p. 559) reads: 

F(a, b; c; z) = (1 - z)~ a F (a, c - 6; c; -^j . (14) 

Using equation (fli"]) for the hypergeometric series of U\{r]) we get: 

2 Fi (l,-iS;2-iS--e- 21 ^) = 2 Fi l,2;2-iS; 

v ' ' ' 2 cos 77/ V 1 + e 1V 

Thus: 

Ui(v) = 9 g, 1 - f 1, 2; 2 - iS; 1 ) . (15) 
2e ,: "? cos 777 \ 1 + / 

On the other hand, we can use the following relationship: 

(l-z)-P = 2 F 1 (p,b;b;z) 
for p = 1 and S = in f)15[) . which leads to: 

Wi(n) = 2 cos Tfae-*^ = e"^, 

2 cos 77/ 

as expected. 
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